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Key diagnostics of non-thermal plasmas, such as the position of the syn- 
chrotron self-absorption turnover or the Faraday rotation optical depth, depend 
crucially on the unobserved low-energy electron spectrum. We investigate the 
^ ' effect of physically viable assumptions about the low-energy electron spectrum 

on the magnetic field strength as inferred from the position of the peak in the 
' synchrotron spectrum marking the transition from optically thick to optically 

. thin radiation. To this end the kinetic equations are written down and sta- 

' tionary solutions are given. In this context the inviability of the delta function 

approximation of the synchrotron power spectrum is discussed. Using two dif- 
ferent approaches, we show that the field determination is independent of the 
assumptions about the low energy electron spectrum within an order of mag- 
On' nitude. As an example, we use our result to obtain a small correction to the 
q limiting brightness temperature. 

1 Introduction 

The evaluation of the magnetic field strength of radio sources may in principle 
■ be based on the fact that their synchrotron spectra exhibit a maximum, which 

arises due to self absorption in the emitting plasma (Rybicki & Lightman 1979). 
This necessitates the full knowledge of the electron^] distribution. However, 
whereas the distribution in the optically thin regime can be inferred from the 
spectrum and is known to be of power law form, it remains uncertain in the 
optically thick regime, so that for the latter some assumption is needed. Obvi- 
ously, the simplest alternative would be to assume an unbroken power law (cf. 
Pacholczyk 1970). 

In this article, the differences between an unbroken power law and some 
other physically more realistic distributions are investigated. Throughout, apart 
from the assumption of homogeneity and isotropy, the discussion will be kept 
as independent as possible from specific source parameters. In Section |^ the 
kinetic equations governing emission and absorption in a synchrotron plasma 

1 Although not stated explicitly, all results of this paper are valid for electron and electron- 
positron plasmas alike. 
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are reviewed, and two simple steady state solutions are given. They motivate 
the use of either a power law with index -3 or a thermal electron distribution 
in the optically thick regime. Section || deals with the magnetic field strengths 
one obtains under the assumption that the transition from the optically thick 
to the optically thin regime corresponds to the turnover energy in the electron 
distribution. Piecewise power laws and thermal distributions for low energies 
are discussed. In Section |^ the determination of the magnetic field strength 
from the maximum intensity of the spectrum is investigated. The uncertainty 
in the limiting brightness temperature is discussed briefly in Section [| 



2 Kinetic equations 

Treating synchrotron radiation as a special case of bremsstrahlung, we can see 
that synchrotron emission as a spontaneous process must be accompanied by the 
corresponding induced processes. Accordingly, it may be described by means 
of the Einstein coefficients, so that the Einstein relations A21 — ^%-B>2i and 
B12 = B21 may be applied. In order to simplify, in the following we assume 
the electron energies e are highly relativistic with e Huj (w being the angu- 
lar frequency of an emitted photon) and that both the magnetic field and the 
electron distribution may be treated as isotropic. Replacing A 2 i by the power 
spectrum P(e, o>) and employing the Einstein relations then yields the kinetic 
equations describing a synchrotron emitting plasma (McCray 1969, Norman & 
ter Haar 1974): 

8W{u) /■.,-=. . f , tt 2 c 3 2 8 //(e) 



dt 



In these formulae, / stands for the electron density per energy e, W for the 
photon energy density per angular frequency u>. P is the angular average of P, 
i.e. P = J i7T Pdfl. The kinetic equations obtained in this manner remain valid 
in case of a homogeneous magnetic field and an isotropic photon distribution. In 
the stationary case comparison of the kinetic equations and the transfer equation 
^-j- = j v —a^v with emission coefficient j„, absorption coefficient a v , and path 
length s yields an expression for ot v : 



Turning now from these rather general considerations to the specific exam- 
ple of synchrotron radiation, one starts by defining x = uj/lu c and F(x) = 
x J°° K 5 / 3 (£_)d£ l , where uj c = 3e 2 ePsina/2mj?c 5 with the elementary charge 
e, the electron mass m e and the pitch angle a between magnetic field and 
electron velocity. K 5 / 3 constitutes the modified Bessel function of order 5/3. 
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Using these definitions, one may express the synchrotron power spectrum as 
P(e,u) = {^e z B sina/2nm 1 c 2 ) 1 l 2 F{x) (Ginzburg & Syrovatskii 1965). Its 
angular average is given by P(e,cj) = Q(x) P (e, w)| sino=1 with x = x\ sina=1 , 
where for 

the approximation Q(x) w £^ =0 Q„£ n = 0.822533-0.181121 x+0.0370118 i 2 - 

0. 003166231 3 can be used. (Because of P(e,u) < P(e, w)| sinQ , =1 , < 1, so 
that for i — ► oo the approximation becomes arbitrarily inaccurate. However, 
Q{x) will always be part of an integrand containing e~ x , which renders this inac- 
curacy irrelevant.) Concerning stationary distributions, two important solutions 
of the kinetic equations can be stated (Norman & ter Haar 1974). Firstly, the 
distributions may be those of thermal equilibrium in the Rayleigh- Jeans limit, 

1. e. /(e) = /oexp (— e/7r 2 c 3 f2o), W(ui) = O w 2 , where </>o and f2o are arbitrary 
positive constants. Secondly, the only viable power law for the electron dis- 
tribution has the index -3, so that /(e) = /o e 3 - The corresponding photon 

distribution has the form W{lo) = 572 y^^^ 5 ^ 2 - Finally, one may note that 
the form of the kinetic equations precludes the use of the delta approximation 
P(e,oj) = Ci€ 2 S(uj — C2e 2 ) (C\, C2 being constants) for the power spectrum, as, 
despite of yielding the correct photon distributions, this approximation allows 
any electron distribution, which contradicts the fact that only one power law 
index is allowed. 



3 Magnetic field strength 

A typical synchrotron spectrum contains a turnover, marking the transition 
from optically thick to optically thin emission. Whereas, from observation, the 
optically thin emission (and thus the electron spectrum for sufficiently high 
energies) can generally be taken to be a power law with some index P2, the 
precise form of the optically thick spectrum is hard to detect. Therefore the 
electron spectrum is taken to be of the form 



m = 



F(e) , e < e T 
/oe P2 , £>£T 



where various functions F will be discussed. The simplest choice is to assume 
an unbroken power law spectrum, i.e. F(e) = /oe P2 . Quantities relating to 
this choice are denoted by the index a. Inserting / a (e) and the (averaged) 
synchrotron power spectrum P(e, u>) into equation (Q), one obtains, using the 
relation / x^F(x)dx = (f + |) T (f + §) with fj, > -f (Westfold 1959), 



3 

rtnfnR. 2 " — " 2 

n=0 



ao/oPa 2 v 2 2 /i(p 2 ,n), (2) 
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where a = R = -t^tt and 

u S7rm c 5 47rmgc 5 

h( P , n) = Q n [l- |) n + 1 _| r (6) r (6) » 

where the & are defined as £i = (6n + 22 — 3p)/12 and £2 = (6n + 2 — 3p)/12. 
Note that equation can easily be solved for the magnetic field strength B. A 
second, slightly more general choice for /, which will be indicated by the index b, 
is a piecewise power law, so that F(e) = /beS? ~ Pl e Pl . Here the constant dp 2-Pl 
is implied by demanding continuity in ex. Using the approximation F(x) ~ 
Foi^e-* = 1.8i 3 e- 5 (Melrose 1980) together with x{e,v) = x\ s - na=l = v/Re 2 
and 

g{p,n,x u x 2 ) ee -F Q n {C 3 r| (Ci) - rg (( 2 )} 

(where £1 ee n + cp — p/2, £2 ee n + cf> — p/2+1 and £3 = n + — 1) one gets for 
the ratio ct Ult t,/a U2 , a of the absorption coefficients computed with / a and /b 

Q^i.b _ S»=0 {"dcn(^l, eT,b,Pi,P2, n)} 

where 

adcn(^l,eT,b,Pl,P2,«) 

_ P2_ P2 _2 

= ao,b/o^? b 2 ^l 2 .9(^2, n,0,x(^i,e T ,b)) 

PI PI o 

+ao,b/o£ P2_Pl -R b 2 ^l 2 .9bi;"' J S( !/ i:eT,b),oo). 

Now the synchrotron power spectrum takes its maximum at the frequency 
u mp (e) = 0.29fi?e 2 , so that i(^ mp (e T ), e T ) = 0.29f i?e|/i?4 = 0.29f. As 
ex constitutes the transition between the parts of the electron spectrum yield- 
ing optically thin and optically thick radiation, it is reasonable to assume that 
the optical depth corresponding to i^mp(eT) is approximately 1, i.e. that for a 
homogeneous plasma a„ mp ( £T )Z/ w 1, where L denotes the spatial extension of 
the plasma. (In the case of an unbroken power law, this condition is used to 
define et-) Hence 1 ss Q ? i/ mp (e T ),b/< :k y mp (cT),a- If ° ne allows for different turnover 
energies, this may be generalized to 1 w a„ mp ( eT b ),b/ Q! i/ mp (£T a ).a> an d using 
equation @ and i(t / mp (eT), ct) == 0.29^ together with the ratio r b = £T,b/eT.a 
one can obtain the ratio 6b = B\> /-B a of the magnetic field strengths for / a and 
/ b , respectively: 

, P2 4 ELo{7 1 + (^) i(Pl - P2) 7 2 } 
»b = r b ' 3 

Here the ji denote ji — g(p 2 ,n,0,^^-Tr) and 72 = g(pi,n,^^w,oc). If the 
turnover frequency is independent of the electron spectrum (i.e. if i / m p,a(eT,a) — 
i/ mp b ( e T,b))j the relation br 2 = 1 holds valid. Hence we assume that br 2 = N 
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Figure 1 : Ratio 6b of magnetic field strengths due to /b and / a as a function of 
the power law index p2 for N — 1 (solid), 0.5 (long dashed) and 2 (dot-dashed 
line) . 



with some arbitrary, but constant N. In view of the discussion at the end of 
the previous section, one should take p\ = —3. For 7V=0.25, 1, and 2, the plots 
of 6b versus p2 corresponding to this choice of p\ are shown in Fig. [l| Note 
that 6b(p2 = —3) 7^ 1 for N = 1 results from the use of the approximation 
F(x) « 1.8x 03 e~ x in computing a Vj b- A simple generalization of /b is achieved 
by dropping the assumption that the electron spectrum extends down to zero 
energy (which, due to the finite electron rest energy, is questionable anyway), 
i.e. by taking a spectrum /b'(e) = /b(e)@(e— eo) rather than /b(e), where eo > 
and denotes the Heaviside function. This merely involves replacing x-i = oo 
in the argument of g(p, n, x%,X2 = oo) by x(v, eo) with the cutoff energy eo. If 
one introduces the ratio e = ex.b/eoi one obtains x(v mp ( 
and thus for the ratio 6b' = / £? a 



,(e T ,b),eo) = 0.29fe 2 



„P2-4 



i f 0.297T \ : k(Pi-P2) - \ 
2^n=0 \7l + {— —) 72 1 



where 72 = g(pi,n, ^|^7r, 0.29f e 2 ) and 71 is defined as above. Again, pi should 
be set to -3. The plots of 6b' as a function of e resulting from this formula 
are given in Fig. for several values of P2 and N = 1. As expected, in the 
limit of e — > 00 (no energy cutoff) 6b' approaches 6b; for values of e exceed- 
ing 10, the two are essentially indistinguishable. Referring to the examples 
in the previous section, it is tempting to assume, as a further possibility for 
the electron spectrum, a thermal distribution for low energies, i.e. to choose 
F(e) = / e /5£T e^ 2 " 2 e 2 e- /3£ , where /3 = l/k B T with the Boltzmann constant fc B 
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Figure 2: Ratio of magnetic field strengths due to /t/ and / a as a function 
of the ratio e of turnover and cutoff energy of fb> ■ N — 1 is assumed, and the 
lines are given for the power law indexes p 2 = — 5 (solid), -3 (dot-dashed), -2 
(long dashed), and (dotted line). 



and the temperature T. One can see easily that continuity in £t is ensured. The 
quantities relating to this choice of /(e) are marked by the index c. Inserting 
/ c (e) into equation ([I]) and using i(i/ mp (ex,c), £t,c) = 0.29^ yields 



3 

£ 

n=0 



|«0,c/o 



2 T-.-g-—2, 



7i 



where A(pi,p 2 ) 



+ a ,JoA(p 1 ,p 2 )f3e^^g(f3e T , c ,n)} 

P2 i i 

'0.29ir\ - — + 1 



R c 2 2 
dxx 



V 2 2 7/ 



^m P (eT,c), and 



Now in order to simplify the discussion, we assume Nth to be minimized. This 
assumption is motivated by the low photon density in the optically thin part of 
jet spectra, as implied by polarisation measurements (Jones 1988). Hence one 
0, which together with u = /3ct, c leads to 



demands 4^ 



- 2(p 2 - 2)(e" - l)u- A + 2(e" - (p 2 - l))u^ - p 2 u~\ 

Let u m i n = u m in(p2) be the solution of this implicit equation. Then using Eqs. 
and and defining 7 = g(u m i n , n), one obtains for the ratio b c = B c /B & 



ELo{7i + (»)^ +i Um inC-» 7} 



J2n=0 h (P2, 
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Figure 3: Ratio b c of magnetic field strengths due to f c and / a as a function of 
the power law index p2 for N — 1 (solid), 0.5 (long dashed) and 2 (dot-dashed 
line) . 

where r c = £T,c/eT,a, the analogue of rt,, is assumed to be given by br\ = N with 
some constant N. The plots of b c versus P2 for several values of N are shown 
in Fig. H Note that the restriction to values of P2 less than -1 is necessitated 
by the demand that u m i n be positive. 

4 Maximum intensity 

So far, it has been assumed that the optical depth of the frequency correspond- 
ing to the turnover energy is equal to 1. However, from an observational point 
of view, it is clearly more desirable to deduce the magnetic field from the posi- 
tion of the maximum intensity. In order to achieve this aim, one assumes the 
radiative transfer to be one-dimensional and starts from the general formula for 
the intensity l v of a homogeneous source (Pacholczyk 1970), 

/„ = ^(l-e-^ L ). (4) 
a v 

Now for a piecewise power law with a lower cutoff energy e~ 1 tT, using the nota- 
tions introduced in the previous section, one may write the absorption coefficient 
as a v = ^2 i=1 ^2 n - aoA i R~^i/'i'~ 2 g{p i ,n,xu,X2i), where in = x(v, e _1 e T ), 
X\2 = x 2 \ = x(v,e T ), x 2 2 = x(v,oo) = 0, A\ = foe^~ Pl , and A 2 = f a . If one 
introduces the turnover frequency vq = !/ mp (eT), the Xik become functions of 
this frequency only, i.e. they may be considered to be independent of the mag- 
netic field. Accordingly, a v takes the form a v , Ma — ^i=i foFai(v, vo)G a i{B) 
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with polynomial functions G a i(B) = B Ki . However, from ao oc R oc B and 
Ai oc e^ 2 ~ Pl oc vl /2 B^-^/ 2 it follows that k x = K 2 = 1 - = k, yielding 

= G a {v,vo)B K and hence ^;Q!i/ = foFQ a (v,VQ)B K . Similarly one obtains 
> = f F j (u,u )B x and |;> = f F d] {v,v )B x (with A = ±(1 - p 2 )) for the 
emissivity and its derivative. 

The intensity takes its maximum when its derivative vanishes. Therefore, 
by differentiating equation (jij) one obtains the condition 

= C 1 K ax ,^ )(l-e- C2 ^— 

+ C 3 ((*w^o))ze- c ^— ■ ^, (5) 

where x = / £B K , Ci = F dj /F a -{Fj/F*)F da , C 2 = F a , andC 3 = (Fj/F a )Fg a . 
Note that this equation is formally independent of both the source extension 
and the overall electron density (i.e. of /o). Only when computing the magnetic 
field strength by means of B — (x/L) 1 ^ these parameters reenter. Clearly, 
for two models with the same values of L/o, i-Wx, and vq, the ratio of the 
magnetic field strengths to the power of k equals that of the x obtained from 
equation (||). Hence, it is straightforward to use the ratios 6b and 6b' as de- 
fined in the previous section. Incidentally, this leads to a further simplifica- 
tion: Due to the fact that x{v, e _1 £T) is a function of y = v/vq and that 

A\ oc e^ 2 Pl ' oc v^ V2 , the functions G defined above may be written as 

£2.-2 - ¥2.— 3 - £2 + 1 - 

G a (v,v ) = v 2 G a (y), Gqo.{v,vq) = f 2 G da (y), Gj{v,v a ) = v 2 2 G 3 (y), 

E2. — 1 - 2. - 

and Gdj{v, vq) = v 2 2 Gdj(y), which implies C\(v, v ) = v§C\(y), C 2 {v, v ) = 

P2._2 - E2. — 1. - 

Vq C 2 (y), and 63(^,^0) = v Q 2 2 C(y). Defining x(y) as the solution of 

Ci(y)(l - e -° 2(y)s{y) ) + C 3 {y)x(y)e- C2iy ^ = 

3/2 

and multiplying both sides of this equation by v , one can see easily that 

i-j _£2._|_2_ 

equation (g) is solved by x(^ max , v a ) = !/ 2 x(y max ), where y max = t'max/^o- 
Accordingly, the ratios :rb / £a and Xb' / x a of the x computed for the electron 
distributions a, b and b' of the previous section depend on 2/ max only. How- 
ever, (assuming L and /o are the same for all distributions) these ratios clearly 
constitute the ratios of the respective magnetic field strengths to the power of 
K, so that 6 b = 6 b (y max ) = {x h /x a ) 1/K , b h , = b h >(y max ) = {xw/x a ) 1/K . Hence 
the magnetic field ratios can be computed by means of solving equation (J5J) for 
the distributions involved. For several values of p 2 and e the results obtained 
numerically in this way are plotted in Figs. || to ^. 

Now although this choice is not motivated by the discussion of Section |], 
it is instructive to examine an electron distribution with a positive power law 
index p\ in the low energy regime. For the sake of definiteness, pi = 2 is taken 
as an example, which may be regarded as a crude approximation of the thermal 
distribution. One then obtains the ratios 6b plotted in Fig. [?]. Clearly, the result 
resembles that for p\ = — 3 and e = 2, which implies that the contribution of 
the low energy part of the electron distribution is rather small. 
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Figure 4: The ratio 6b as a function of I'max/t'o for several power law indexes p2- 
The lines correspond to P2 = —5 (solid), -3 (dot-dashed), and -2 (long dashed 
line) . 
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Figure 5: The ratio 6b' for e = 4 as a function of u max /vo for the same power 
law indexes as in fig 0. 
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x a and thus B a do not depend on v (see equation @ below). Hence, as- 
suming some specific value for u max , the ratios bi are given by bi(u max /Po) = 
CBi(v max _/vo), where the constant C is independent of the spectrum i. There- 
fore Bi(z\) I Bj(z-i) = bi(zx)/bj(z2) (where z = v max /vo), which allows a direct 
comparison of Figs. f|to[7| This clearly shows that the magnetic fields obtained 
for the various electron spectra discussed in these figures are of the same order 
of magnitude.^ 

Obviously, in order to get Xb (or Xb>) instead of the ratio Xb/x a (or Xb>/x a , 
respectively), one needs to know the value of x a . Now, as is evident from 

_22.-)_2_ 

the above discussion, x a must be of the form x a (v max , vq) = Vq 2 s;(j/max)- 
However, concerning an unbroken power law, the turnover frequency has no 
physical meaning, as there is no turnover at all. Still, it may be retained as a 
formal parameter, so that the distribution is viewed as a piecewise power law, 
where the power happens to be the same for all pieces. This implies that for 
u any value may be assumed. Choosing — u max , one obtains y max = 1 and 
thus 

pi | o 

x a = x(l)v max . (6) 

Finally, an important technical remark should be made: The aim of the foregoing 
discussion was to gain results depending on as few parameters as possible. In 
order to achieve this, the turnover frequency rather than energy was employed. 
However, in a given model, the latter rather than the former will be given. 
In this case the magnetic field ratio cannot simply be read off, and the field 
strength has to be inferred from demanding self-consistency: Assuming some 
specific values of L, /o and the magnetic field B, one can compute the turnover 
frequency from the corresponding energy and thus obtain the ratio x\,/x a (or 
Xb'Jx a ). Together with the value of x a taken from equation (||) one gets the 
value of Xb (or Xb>), leading to the magnetic field strength, which obviously 
should equal the one originally assumed. 



5 An example: peak brightness temperature 

As a straightforward example of using the results obtained so far, we discuss the 
maximum brightness temperature of a (synchrotron) source with an isotropic 
magnetic field (cf. Kellermann 1974). To do so, one starts by noting that the 
ratio of the synchrotron and (to first order) inverse Compton losses is given by 

P c _ "rad _ L/(ttp 2 c) _ 32ir£ gcom R 2 F max v c 

Psyn U B B 2 /(8lT) ~ p 2 cB 2 

Here, w ra d and u b denote the energy density of the radiation and the magnetic 
field, respectively, p the radius of the source, R its distance from the observer, 

2 Obviously, this is not true, if, in Figs. |^ to ^, i^ ma x/^o happens to be smaller than the 
value corresponding to the maximum of the relevant curve. In that case, however, due to the 
steepness of the curve, one should distrust the magnetic field determination anyway. 
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and v c the upper cut-off frequency. Clearly, p — R9/2 with the angular size 
9. The luminosity L is approximated as L rs CgcomAn R 2 F m iy c , where F m is the 
observed maximum flux density and £ gC om takes into account any deviation from 
isotropic radiation. 

The solid angle of the source (as seen from the observer) is given by ir9 2 /A. 
Hence, assuming that the one-dimensional transfer equation (equation (Q)) may 
be applied, one can show from the formulae for the absorption and emission 
coefficient that for an unbroken power law electron distribution with index p 
the magnetic field strength may be expressed as 



B pl = 2.488 x 10 - 62 ( ^4 ) ! I 



</(^ 2 ) 



1 erg /cm 



(^max \ 1 
1 Hz/ 



G. 



T max constitutes the optical density corresponding to the frequency ^ max of the 
maximum intensity. The functions a(jp) and b{p) are of the order 1 and can be 
found in Longair (1994). The magnetic field actually present can be written as 



B = bB 



pi, 



where the factor b is dependent on the electron spectrum (cf. Sections ^ and | 
In addition, the peak brightness temperature of the source is defined as 



c 2 F 

Ti L - ± max 

max 



2fc B i& aac tt^/4- 

We now investigate the influence of the low energy electron spectrum on the 
limiting brightness temperature T^ ax ~ 10 12 K for the onset of the inverse 
Compton catastrophe, i.e. the temperature T max for which P c /P S y n = 1. Com- 
bining the formulae obtained so far, including second order inverse Compton 
losses, and choosing the power law index p — —3 yields 

n r r» Cgcom ( ^ma) 



= 0.52- 



P syn ' b 2 \10 12 KJ V 1 MHz 



n. rp Sgcom / 1 max \ / v c \ 

b 2 ll0 12 K/ VI MHz/ 



If we assume that 0.3 < b = B/B p \ < 3, we see that the values of Tj£ 
sponding to B v \ and B differ by a factor less than approximately 3 2 / 5 w 1.6. 
Hence, the uncertainty in the magnetic field determination due to the unknown 
optically thick electron spectrum has a rather small impact on the limiting 
brightness temperature. Indeed, the uncertainty in the source geometry (i.e. in 
Cgoom) may be considerably more important. 
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6 Conclusions 



Based on the discussion of stationary solutions to the kinetic equations in a 
homogeneous and isotropic synchrotron plasma, the influence of the optically 
thick part of the electron spectrum on the magnetic field determination from 
the position of the turnover in the photon spectrum has been investigated. By 
employing the assumption that the optical depth equals unity at the turnover 
frequency it was shown that the error made by choosing some specific form 
for the low energy electron spectrum can be estimated to be within one order 
of magnitude. For power law distributions this result was corroborated by a 
more rigorous treatment. All results obtained are independent from the spatial 
extension of the source and the overall number of electrons. They show, for 
example, that the limiting brightness temperature is influenced to a small extent 
only. 

Thus we have seen that that determining a magnetic field strength from the 
turnover of a synchrotron spectrum yields a value that is of the correct order 
of magnitude. Of course, in actual computations other uncertainties such as in 
the source geometry must also be taken into account. 
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